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We define the notion of Bernstein measures and Bernstein approximations over 
general convex polytopes. This generahzes well-known Bernstein polynomials which 
are used to prove the Weierstrass approximation theorem on one dimensional intervals. 
We discuss some properties of Bernstein measures and approximations, and prove an 
asymptotic expansion of the Bernstein approximations for smooth functions which 
is a generalization of the asymptotic expansion of the Bernstein polynomials on the 
standard m-simplex obtained by Abel-Ivan and Hormander. These are different from 
the Bergman-Bernstein approximations over Delzant polytopes recently introduced by 
Zelditch. We discuss relations between Bernstein approximations defined in this paper 
and Zelditch's Bergman-Bernstein approximations. 



1 Introduction 

It is quite well-known that Bernstein ([B]) introduced the polynomial 



to prove the Weierstrass approximation theorem. After his work, several properties and 
generalizations of the Bernstein polynomials ([1]) are considered. For example, the Bernstein 
polynomials on the standard m-dimensional simplex P is defined by 



where / is a function on P, a = (ai, . . . , am) G Z>o and ||a|| = X^jli ^j- The reader may 
be referred to the book [L] in which one can find many interesting properties and some 
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generalizations of the Bernstein polynomials. Some relations between Bernstein polynomials 
and Brownian motions are discussed in [Kj. 

Rather recently, Abel-Ivan ( |AI] ) and Hormander ( |Ho2] ) obtained a complete asymptotic 
expansion of the Bernstein polynomials on the standard m-dimensional simplex defined by 
([2]). In particular, Hormander used it to analyze asymptotic behavior of certain Bergman 
kernels near the boundary of the unit ball in a complex space. 

Further generalization was considered by Zelditch (|Z|, see also ^]), who defined the 
Bergman-Bernstein approximations on regular simple lattice {i.e., Delzant) polytopes. He 
defined it by using Bergman-Szego kernels on a toric Kahler manifold whose moment poly- 
tope is a given Delzant polytope. He obtained an asymptotic expansion of the Bergman- 
Bernstein approximations by using the theory of Toeplitz operators and the harmonic analy- 
sis on toric Kahler manifolds developed in {STZlj , |STZ2] , [SoZlj , and applied the Bergman- 
Bernstein approximations to obtain an asymptotic expansion of the Dedekind-Riemann sum 



over the Delzant polytope P in M*". Note that an asymptotic expansion of the Dedekind- 
Riemann sum ([3]), for simple lattice polytopes, in a form similar to formulas of Euler- 
Maclaurin type was obtained by Guillemin- Sternberg [iG'Sj . 

In view of this, it would be interesting to consider some generalizations of the Bernstein 
polynomials for general convex polytopes and its applications to various directions. Our aim 
in this paper is to define in a natural way a generalization of the Bernstein polynomials on 
general convex polytopes and obtain its asymptotic expansion. We avoid to use harmonic 
analysis on toric varieties to define the Bernstein approximations on general polytopes, be- 
cause the toric variety corresponding to a general polytope will have singularities. Instead, 
we use the setting-up and some analysis discussed in |TZ] . Our main theorems are described 
in Section [2] and proved in Section [31 

Our Bernstein approximations defined in Section [2] are different from Bergman-Bernstein 
approximations defined in [Z]. In [Z], Zelditch used a general Kahler metric on a toric 
Kahler manifold. One may say that, when the polytope is Delzant, we use only the Kahler 
metric induced by the restriction of the Fubini-Study metric through a monomial embed- 
ding. But, our Bernstein approximations can be defined for general polytopes. Recently, 
Song-Zelditch ( |SoZ2] ) shows a large deviations principle for the Bergman-Bernstein mea- 
sures. In particular, they give a concrete description of the rate functions. Our Bernstein 
approximations are defined in terms of dilated convolution powers of Bernstein measures, 
and hence they satisfy a large deviations principle. We give a concrete description of the 
rate functions for finitely supported Bernstein measures in Subsection 13. 4[ which eventually 
coincide with Song-Zelditch rate functions ( |SoZ2] ) for Bergman-Bernstein measures defined 
from Fubini-Study metric on a projective toric manifold when the polytope is Delzant. We 
discuss in Section H] relations between our Bernstein approximations and Bergman-Bernstein 
approximat ions . 

Acknowledgments. The author would like to thank to professor Zelditch for his valu- 
able comments on the topics discussed here, and professor Bando for his pointing out to the 
author on the facts given in Subsection 12.21 
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2 Bernstein measures on convex polytopes 
2.1 Definitions and main theorems 

Let P C be a polytope which has non-empty interior, P°. The barycenter, 

b{fi) = / zdfi{z), (4) 

of a given probabihty measure fi on R"* is contained in the convex hull of the support, 
supp(/i), of /i, and hence it defines a continuous surjection 

b:M{P)^P, (5) 

where Ai{P) denotes the set of probability measures on the polytope P with the weak-* 
topology. We denote C{P) the space of all continuous functions on P. 

Definition 2.1 A continuous section B : P ^ A4{P) of the barycenter map satisfying 
the following two conditions is called the Bernstein measure on the polytope P: 

(1) For any f G C{P), the function B{f), defined by 

B{f){x):= jj{z)dB^{z), dB,:=B{x)eMiP), x E P, (6) 
IS m C°°(P°) nC(P). 

(2) There exists a smooth function K : P" ^ Sym(m, M) such that for any f G C{P), we 
have 

VB{f){x) = f{z)K{x){z - x) dB^iz), X G P^ 

where Sym(m, R) denotes the space of all symmetric m xm real matrices. We call the 
function K : P° ^ Sym(m, M) the defining matrix of the Bernstein measure B. 

Definition 2.2 Let B : P ^ M{P) be a Bernstein measure on P. For any f G C(P) and 
any positive integer N, we define the function B^i^f) on P by 

B^{f){x):= Jj{z)dB^{z), (7) 

where dB^ is the probability measure on P given by 

dB^ = {DyN)*{dB,*---*dB^), Di/^ : R'" 9 X ^ x/AT G R'". (8) 

We call the function B^if) the Bernstein approximation of f . 

Concrete expression for the dilated convolution measure dB^ of B{x) defined in ^ is 
given by 

/ f{z)dB^A^)= [ f{{z^'^ + --- + z^'^^)/N)dB,iz^'^)---dB,iz^''^), (9) 

Jp Jpx-xP 
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where the integral domain P x ■ • • x P is the times product of the polytope P. 

To our knowledge, generalizations of the classical Bernstein polynomials previously con- 
sidered are defined by using finitely supported probability measures. This is true also for the 
Bergman-Bernstein approximation in [Z]. We do not assume here that the measure dB^ to 
have finite support. However, a Bernstein measure having a finite support in the following 
sense might be important. 

Definition 2.3 Let B : P ^ Ai{P) be a Bernstein measure. Then, B is called a finitely 
supported Bernstein measure if there exists a finite set S in P such that supp(P(a;)) = S for 
any x G P° and supp(P(x)) C S for any x G P. The finite set S is called the support of B. 

Now, let ;B : P ^ Ai{P) be a finitely supported Bernstein measure with support S G P. 
Then, it is easy to show that, for each a E S, there is a function nia G C°°(P°) fl C(P) such 
that dBx = ^Q,g5 "T'a(3^)^a for each x G P. The functions rUa, a E S satisfy the following: 

(1) rria is positive on P°, and we have J2aes "^a(a^) = 1 for all x G P; 

(2) for each x G P, we have '^a&s^a{x)a = x, which implies that the convex hull of S is 



(3) we have VmQ,(x) = ma{x)K{x){a — x) for each x G P". 

Conversely, a finite set 5* C P whose convex hull is P and functions rria G C°°{P°) fl C(P) 
for each a E S satisfying the above three conditions define a Bernstein measure B{x) = 

The following lemma gives a concrete expression for the dilated convolution powers dB^ 
for a given finitely supported Bernstein measure dBx- 

Lemma 2.4 Let B : P —>■ Ai{P) be a finitely supported Bernstein measure with the support 
S, and write it as B{x) = J2aes^»(-'')^»- Then, the probability measure dB^ is written as 




(10) 



where Sn C NP is the finite set defined by 



SN = {ieNP;j 



Pi-l \- Pn for some Pi, . . . , Pn G S} 



(11) 



and the function m]^ G C (P) fl C' 



'oo 



(P°) is defined by 





(12) 



/3i,...,/3jvGS 

/3i+---+/3iv=7 



Proof. By the Bernstein approximation B^if) is given by 



BNm^) = Yl ^((^1 + ■ ■ ■ + PN)/N)mp,{x) ■ ■ ■m^^(x) 



/3i,...,/3jvG5 



= /(7/A^)"^/3i(a;)---m^^(x), 



TeSjv /3l,..../3iV 

/3i+---+/3iv=7 
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which shows the equations (ITOl) . (fT2|) . □ 

We mention some remarks on the definition of Bernstein measures and Bernstein approx- 
imations. 

Remark 2.5 It is easy to show that the Bernstein approximation B]sr{f) for / G C{P) 
with respect to a Bernstein measure B . P ^ M.{P) is in C°°{P°) (1 C{P). The condition 
(2) in Definition 12.11 comes from Hormander's proof ( |Ho2] ) of an asymptotic expansion of 
the Bernstein polynomials ([2]) on an m-dimensional simplex. In fact, he uses a differential- 
recurrence formula for a family of functions (which is, in our setting-up, the function I^^a 
defined in (l29l)). In his proof of the recurrence formula, he uses the matrix 'A(x)' given 
in Example 12.61 below. The recurrence formula for the functions I^^a defined in (!29|) is 
given in Lemma 13. 7[ The condition (2) in Definition 12.11 assures that the functions lN,a 
are polynomials in (Lemma 13.81) . and it provides a computable representation for each 
differential operator in the asymptotic expansion (1211) in Theorem 12. 9[ 

In general, for a given section B : P ^ Ai{P) of the barycenter map define the 
probability measure dB^ by dH]). Then by the Law of Large Numbers, dB^ tends weakly 
to the Dirac measure 5^; at x G P. Furthermore, this convergence is uniform on P. See 
Subsection 12.21 We also note that, since the measures dB^ is defined as a dilated convolution 
powers, it satisfies, for example, the Central Limit Theorem (with a suitable dilation) and a 
large deviations principle at least for fixed x. In particular, the rate functions for the large 
deviations principle for the measure dB!^ are given in Subsection 13.41 Proposition 13.161 

The measure dB!^ is also a section of the barycenter map, that is the barycenter of 
dB^ is X. But, it might not be necessary to define dB^ as a dilated convolution power. 
One might need only the properties that the measure dB^ has some regularity in x G 
and its barycenter suitably converges to x. In fact, the measure defining the Bergman- 
Bernstein approximation in ^ is not a section of the barycenter map, but it has these 
properties. However, as in Section [3l that the measure dB^ is a dilated convolution power 
makes analysis much easier than something which satisfies only the above properties. This 
is one of the main differences between the Bernstein approximation in this paper and the 
Bergman-Bernstein approximation. See Section H] for details. 

Next, we give some examples of Bernstein measures. 

Example 2.6 Let P be an m-dimensional standard simplex, 

m 

P = {x = (xi, . . .,Xm) e ; Xj > 0, ||x|| := < 1}. 

i=i 

Let S = {cq := 0, ei, . . . , Cm} = P H Z™ where {ej}Y=i is the standard basis of M™. Define 
the functions rriej G C°°{P) (j = 0, 1, . . . , m) by 

meo(x) = 1 - ||x||, me^-(x) = Xj, j = l,...,m. 

Then, it is easy to see that the measure B{x) = X]j=o "^ej(^)'^ej defines a Bernstein measure 
with the Bernstein approximation given by ([2]). The defining matrix K{x), x G P°, is given 
by(m) 

\x, l-||x||;.^. 
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which is the inverse of the matrix 

A(x) = {xjSij — XiXj)-- . 

Example 2.7 // P and Q are polytopes in with P" ^ ^, Q° ^ 0, then, their product 
PxQ is a polytope in R^™. Let dB^ and dB^ he Bernstein measures on P and Q, respectively. 
Then, the product measure dB(^xy) := dB^dB^ is a Bernstein measure on P x Q. For 
example, let B^ : [0, 1] A^([0, 1]) be the Bernstein measure given in the previous example 
(with m = 1). Then the Bernstein approximation associated to the Bernstein measure 

■.= dBl---dBl^, {xu...,xm) e [0,ir 

on [0, 1]™ is given by 

5 Xm ) 

N 



J2 f{h/N,...,kjN)l[(\x^^{l-x.) 

l,...,km=0 1=1 ^ *^ 



N-k, 



This is a well-known generalization of the original Bernstein measure (m). 

There is an example of a Bernstein measure which has a smooth density. We give such 
an example only on the unit interval [0, 1] as a proposition. We define a function /i = /i(r) 
on M by 

T I — e ^ T 

It is easy to show that the function yU in (|T3l) defines a diffeomorphism between M and (0, 1). 
Let r = t{x) : (0, 1) ^ R be the inverse function of /i. Define a smooth function 5{x) on 
(0, 1) by 

5{x) := logx(r(x)) - xt{x), x(^) = (e^ - 1)/^- (14) 
Finally, we define a smooth function p{z,x) on [0, 1] x (0, 1) by 

Proposition 2.8 For any x G [0, 1], define the measure dBx by 

dBx{z) := p{z, x) dz (x G (0, 1)), dBi{z) := 5i, dBo{z) := Sq. (16) 

Then the map B : [0, 1] ^ 7V1([0, 1]) defined by B{x) = dB^ is a Bernstein measure on [0, 1]. 

Proof. Let us find a Bernstein measure dB^ on the unit interval [0, 1] which has the 
form 

dBx{z) = p{z,x)dz, 2G[0, 1], xG(0,1), (17) 

where p{z,x) is a smooth on [0, 1] x (0, 1). If the measure of the form (fT7I) with a smooth 
function p is a Bernstein measure, then the density function p must satisfy the following: 

1 .1 
p{z,x)dz = l, / zp{z,x)dz = x, dxlog p{z, x) = K{x){z — x), (18) 



where K{x) is a positive smooth function on (0, 1) (the defining matrix of B) which we need 
to specify. We take a potential function ip of K on (0, 1), that is ip" = K. Inserting K = ip" 
for the third equation in (fT8|) and solving it as a differential equation for p, we have 

p{z,x) = C(2)e^(^)+^'(^)(^-^) (19) 

with a positive function C{z). For simplicity, we assume that C{z) is constant, C{z) = C. 
Substituting this for the first and second equation in f|T8l) . we have 

^/ _ ^g-V+v(e^' _ 1)^ a.(<^')2 _ C'e-^'^'+^((^'e^' - e^' + 1). 

Eliminating the term Ce~^'^'^''' from the above, we get 

a;r(e^ - 1) = re^ - + 1, r(x) = v?'(x). (20) 

Note that the defining matrix (function) K must be positive on (0, 1) (see Lemma [3 .ip . and 
hence x t— > t{x) is a diffeomorphism from (0, 1) onto its image. Denote its inverse function 
/i = p{t). From fl20|) . the function p{t) must be given by f|T3|) . The defining function i^'(x) 
must be given by K{x) = l//i'(r(x)) = t'{x), x G (0,1). Therefore, the defining function 
K{x) and its integral t{x) are determined uniquely by the equation f[T5]) . Clearly /i(r) = 
(logx(T))' where the function x is defined in f[T^ . Then, we have 5'(x) = —t{x) and hence 
we can take if = —6. From this choice, we have C = 1 and hence the function p{z, x) must be 
the form (fTSl) . One can easily check that, this function p actually satisfies the equations (fT8|) . 
Note that A{x) := K{x)~^ is continuously extended to [0, 1] with A{0) = A{1) = 0. But, 
the function p{z,x) can not be extended to a; = 0, 1 continuously. However, for x G [0, 1], 
define a probability measure dBx on [0, 1] by (IlGil . For any / G C([0, 1]) and x G [0, 1], we 
set B{f){x) = f{z)dBx{z) which is smooth on (0,1). If / G ^^^([0, 1]), an integration 
by parts shows that B{f){x) — » /(I) when x ^ 1 and B{f){x) /(O) when x 0. 
Since supo<^<i |5(/)(x) - 5(5()(x)| < \\f - g\\c([o,i]), f,9 e C([0, 1]), B{f) is continuous 
on [0,1] for any / G C([0, 1]). Hence dBx is a Bernstein measure on [0,1]. By the above 
discussion, one knows that a Bernstein measure dB^ on [0, 1] having a smooth density is 
uniquely determined, and is given by (ITBi) . (IT^ . □ 

One of our main theorems is the following, which generalizes the asymptotic expansion 
of the Bernstein polynomials on standard simplices given by Abel-Ivan |AIj and Hormander 
jHo2] . 

Theorem 2.9 Let B : P ^ Ai{P) be a Bernstein measure on the polytope P. Then, for 
each non-negative integer u, there exists a differential operator Ly{x,d) of order 1v such 
that, for any f G C°°{P), we have the following asymptotic expansion: 

i?^(/)~$^iV-'^i^,(x,9)/, (21) 

where the expansion holds uniformly on P. Lq{x, d) and Li{x, d) are given by the following: 

Lo(x, d)f = fix), L,{x, d)f = ]p:i{A{x)V^f{x)), (22) 

where V^/ is the Hessian of f and A{x) G Sym(m, M) denotes the inverse of the defining 
matrix K{x) G Sym(m, M). The asymptotic expansion fl?I]) can be differentiated any number 
of times and the resulting expansion holds locally uniformly on P° . 
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The matrix K{x) is not assumed to be non-degenerate. However, it follows that K{x) is 
non-degenerate with the inverse matrix A{x) given by 

A{x) = I {z-x)(^{z-x)dB^{z), xeP". (23) 



jp 

Note that, though the defining matrix K{x) is defined only for x in P°, the matrix A{x) in 
( l23l) is defined and continuous on whole P. We have a computable representation for the 
differential operators L^{x,d). See Section [3] for details. 

Next, to explain a construction of finitely supported Bernstein measures, let us prepare 
some notation. Let S C be a finite set whose convex hull is P. Since P is assumed to 
have non-empty interior, the finite set 5* satisfies the following: 

spaniR{a G 5} = M™. (24) 

Fix a positive function c : S —* M>o on S. Define the map fis,c '■ by 

c(a)e<"'^> 



where {u,v) denotes the usual Euclidean inner product on R'". It is well-known (jF]) that 
IJ,s,c is a diffeomorphism from to P°. Its inverse is denoted by r^c '■ P° For each 

a G S", we define the function ms^c,a by 

Then, clearly, ms^c,a is a positive smooth function on P". Furthermore, it is not hard to 
show that m^co can continuously be extended to the boundary of P so that it defines a 
continuous function on P which is smooth on P° (Lemma 13.101) . Then, since ts,c is the 
inverse map of iis,c, these functions satisfy 

y^^5,c,«(a^)a = X (27) 
for each x G P°, and by the continuity, it holds for each x E P. 

Theorem 2.10 Let S d P and c : S ^ ]R>o be as above. Let ms^c,a, ex E S be the functions 
defined by fl2B]) . Then, the section Bs,c '■ P M.{P) defined by 

Bs,c{^) = ^rns^c,a{^)^» (28) 

is a finitely supported Bernstein measure. Conversely, suppose that there is a finitely sup- 
ported Bernstein measure B with support S . Then, there exists a weight function c : S —>■ ]R>o 
such that B = Bsc- 



We discuss some properties of Bernstein measures and give proofs of Theorems 
I2.1UI in Section [21 As we pointed out in Remark 12.51 the measures dB^ satisfy the large 
deviations principle. In particular, we give a concrete description of the rate functions for 
the large deviations in Subsection 13.41 In Section HI we discuss a difference between the 
Bernstein measures and the Bergman-Bernstein measures defining the Bergman-Bernstein 
approximations introduced in [Z]. In particular, we give some conditions in Proposition 
for when Zelditch's Bergman-Bernstein measures coincide with our Bernstein measures. 
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2.2 General properties of sections of the barycenter map 

Before proceeding to the discussion on Bernstein measures, we give some accounts on prop- 
erties of general continuous sections of the barycenter map b : Ai{P) P on the polytope 
P. Throughout the paper, {ej}^^-^ denotes the standard basis for M*". 

Lemma 2.11 Suppose that we are given a sequence of bounded linear maps Bn '■ C{P) 
C (P) satisfying the following: 

(1) -BAr(l) = 1, B]\f{f) > for non-negative continuous functions f ; 

(2) Bn{x°') — > x°' uniformly as N ^ oo for a G Z>q with \\a\\ < 2. 
Then, Bj\f{f) f uniformly as N ^ oo. 

Proof. First, we note that the following holds: 

• Bjq{f) is real-valued if / is real; 

• |Piv(/)| < 5^(1/1), /eC^(P); 

• \\BNU)\\ciP)<\\f\\ciP),f^C{P). 

Let / G C°°{P). We fix a; G P. For any multi-index a G Z'>q, we set Xa{z) = z"-. Then, 
inserting the Taylor expansion 

f{z) = f{x) + Ri,,{z) + R2,^{z), 

m 

Rl,x{z) = { Vf{x), Z - X) = ^§^{x){Zj - Xj), 

i=i 

i?2,xW= / il-t){V''f{x + tiz-x))iz-x),{z-x))dt 
Jo 

around x G P to B^i^f), we have B^^f) — f = B^^Ri^x) + Bn{R2,x)- By assumption, we 
have Bn{Ri^x){x) = YlJLi '§^{^){BN{Xej){x) — Xj) —>■ uniformly in x G P. Now, we have 
\R2,x{^)\ < \\f\\c\P)\z - and hence 

m 

\Bn{R2,.){x)\ < ||/||c-2(p)5^(Pjv(X2e,)(x)-2x,P^(Xe,)(x)+x2). 

i=i 

Since PAr(x2ej) and B^i^Xej) converge uniformly to x^, Xj, respectively, we conclude that 
BNif) f uniformly on P for any / G C°^{P). Since C°°{P) is dense in C{P), we 
conclude the assertion. □ 

Lemma 2.12 Let B : P —>■ Ai{P) be a continuous section of b : Ai{P) P. For positive 
integer N, we define the measure dB^ by ([8]) and, for f G C{P), the function Pjv(/) G C{P) 
by dZ]). Then, Bj~^{f) converges uniformly to f as N oo for all f G C{P). 
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The following proof using Lemma [2. Ill is pointed out to the author by professor Bando. 

Proof. Note that we have h{dB^) = x. Thus, by Lemma [2. IH we only need to show 
that B^lxa) converges to Xa for each a with ||a|| = 2. Then, by definition of the convolution 
powers (jni), we have 

N 

k=i JP^ - >^P 

Since the second term in the above is 0{1/N), we conclude the assertion. □ 

By using Lemmas 12.111 I2.12[ we further discuss as follows. Let B : P ^ (P) be a 
continuous section of b : Ai{P) P. For a positive integer N, we define the measure dB^ 
by dH]) and for / G C{P) the function B^if) G C{P) by ([7]). For any multi-index a G Z^q, 
we define the function lN,a{x) by 

InA^) ■= / - dB^i^)- (29) 



In the definition (12^ of the functions lN,a, the integral is multiplied by A^H^H. The reason 
for this would be clarified in Section [31 Now, for any / G C°°(P), substituting the Taylor 
expansion 



||a||<2n— 1 ||Q:||=2n 

R^^{z, x) = 2n [ (1 - t)2"-i/W (x + t{z - x)) dt, 



we have 



(30) 



BN{f){x)= E ^[;^^"/(^) + ^2n,iv(x), (31) 
||a||<2n-l 



where the function S2n,N{x) is given by 

S2n,N{^)= E ^ / R2n{z,x){z-XrdB^{z). (32) 



\\a\\=2n 



By Lemma [2.12| ||/Af,Q || c(P) = o(A^""")- Therefore, we have 

\s2nA^)\<\\f\\c^-iP) -J \i^-^r\dB^i^) 

II II n JP 



||a||=2n 



(33) 



which is clearly of order o(l) as N oo uniformly in x G P. However we note that, in the 
expression fl?Il) . each term with > 1 is of order o(l). Thus, one need to find asymptotic 
behavior of the functions /tv,^ to make flSTl) an asymptotic expansion. 
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3 Properties of Bernstein measures; proofs of main 
theorems 



In this section, we discuss some properties of Bernstein measures defined in the previous 
section, and give proofs of main Theorems I2.9[ I2.10[ 

3.1 Properties of Bernstein measures 

Let B . P ^ Ai{P) be a Bernstein measure, and, for any / G C{P) and positive integer 
A^, let -Biv(/) be the Bernstein approximation defined by ([7]). The corresponding defining 
matrix is denoted hy K : P° Sym(m, M). 

Lemma 3.1 For each x G P°, the linear map A{x) : —>■ M™" defined by (1231) is the inverse 
of the symmetric linear map of K{x) : M™. Furthermore, A{x), and hence K{x) is 

positive definite for each x E P° . 

Proof. Differentiating the identity x = JpZdB^iz), we find that, for any u G M*", 



which equals A{x)K{x)u by definition (l23ll of A{x). The fact that A{x) is non-negative (and 
hence positive definite) follows from the definition fl23|) . □ 

Lemma 3.2 For any f e C{P), we have 



J p 

Proof. The assertion is shown by computing VB]s[{f){x) with the definition of the 
dilated convolution power ([9]) in the following way: 



u = 






Px-xP 




fiizi + ■■■ + zm)IN)K{x){z, - x) dB,{zi) ■ ■ ■ dB,{zM) 



X K{x){zi H \- zn 



Nx) dB,{z,) ■ ■ ■ dB^izN) 




*B,){w) 



□ 



Lemma 3.3 We have A{x) = N J {z - 



x)(^{z-x) dB^{z). 
11 



Proof. Since the barycenter of dS!^ is also x, we have x = JpZ dB^ {z). Differentiating 
this identity and using Lemma I3.2[ we obtain 

Id = AT j z® K{x){z - x) dB^{z) = N j {z - x) ® K{x){z - x) dB^{z), 

which means the desired formula. □ 

For any u G M*", we define the first order differential operator Du on P by 

DJ{x) := {A{x)Vf{x),u) = {Vf{x),A{x)u), (34) 

where V/(x) is the gradient of a function / on P. For any smooth function L : P° ^ 
Sym(m, M), we write 



['^L{x)v]u := ^ 



L{x + tv)u, X E P°, u,v & 



t=o 



Lemma 3.4 For any u,v e M™, we have [Du,Dy] = 0. 
Proof. For any A G M™, a simple computation shows 

( VDJix), A ) = ( VV(x)A, A{x)u ) + ( V/(x), [VA{x)\]u ), 
[\/A{x)X]u = I {z — x,u){K{x){z — x), X){z — x) dBx{z), ^ ^ 



where V^/(x) is the Hessian of /. Thus, by setting A = A{x)v, we find 

= {V^fix)A{x)v,A{x)u) + j {z - x,u){z - x,v){Vfix),z - x) dB^{z), 

which is obviously symmetric in u and v. □ 

The following lemma shows an integrability of the defining matrix K, which will be used 
to prove Theorem 12. 10[ 

Lemma 3.5 We have [VK{x)u\v = [VK{x)v]u for any u,v e and x G P°. 
Proof. For x G P°, we have K{x) = A{x)~^, and hence, by fl35|) . we obtain 
[VK{x)u]v = -K{x)[VA{x)u]K{x)v 

= - j {K{x){z - x),v){K{x){z - x),u)K{x){z - x) dB^{z), 

which is clearly symmetric in u and v. □ 

The next proposition shows a uniqueness of the Bernstein measure for a given function 
K : P° ^ Sym(m, M). For any multi-index a = (ai, . . . , «„) G Z>q, we write 

D":=Dr---^r, D,:=D,^, (36) 

where D^^ is defined in ([31]). 
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Proposition 3.6 Let K : P° Sym(m,M"') be a smooth map. Then, the Bernstein 
measure with the defining matrix K is, if it exists, unique. 

Proof. Let B : P M.{P) be a Bernstein measure with the defining matrix K : P° ^ 
Sym(m, R). Then, the matrix K{x) is non-degenerate and the inverse matrix is given by 
A : P ^ Sym(m, M). We note that the inverse matrix A{x) defines the differential operators 
Lj := Xj + Dj, j = 1, . . . ,m, where Xj is the multiphcation operator: Xjf{x) = Xjf{x). 
Then, it is easy to see that 

L,B{f){x)=B{X,f), j = l,...,m, (37) 

where B{f) is defined in ([6]). For each a G Z™q, we denote the monomial with the weight a 
by Xa{z) = z". Then, by f l37|) . we have LjB{xa) = B{Xa+ej)- A direct computation shows 
that 

DkXj - XjDk = { A{x)ej, ) = Ii,ej+e^, (38) 

and which is symmetric in j and k. Hence [Lj,Lk] = 0. We denote = L"^ ■ ■ -L^™ for 
a = (ai, . . . , am) G Z>o- We then obtain 

B{xa) = L"B{1) = ■ 1 (39) 

for each a e Z^g. By ([39]), the restriction of the map B : C{P) C{P) to the set of 
polynomials on P is determined by the matrix K{x) because the differential operators Lj 
are defined in terms of A{x) = K{x)^^. Since B : C{P) C{P) is continuous, it is 
determined by K{x). □ 

Now, for each a G Z>q and positive integer iV, define the function 1^^^ £ C{P) fl C°°(P°) 
by fl29l) with respect to the measure dB^ determined by the Bernstein measure B : P ^ 
M{P). Note that, by definition, 

Ll,ei+ej{x) = { A{x)ei, ^j) = j i^i " " ^j) d^x{z), 

where Zj denotes the j-th coordinate of 2; G M™. From this, we have h^ei+cj G C(P)nC°°(P"). 
Lemma 3.7 The following identities hold; 

(1) /Ar,o(a;) = 1, lN,e,{x) = 0, lN,ei+ej{x) = N { A{x)ei, Cj ) , 

lN,e,+e,+ekix) = N Ii^e.+e.+eki^) for any I <i,j,k <m and N. 

(2) For any a G Z*q, 1 < j < m and N, 

m 

lN,a+ejix) = DjlN,a{x) + ^ aj N ,a-eM)I N ,e,+e,{x) ■ (40) 

i=l 

Proof. First two of (1) are obvious by definition (l29l) . The third formula of (1) follows 
from Lemma [3. 3[ The fourth formula in (1) and (2) follow from direct computations. □ 

By Lemma [3.71 we have D^I^^a G C{P) fl C°°{P°) for any positive integer and multi- 
indices a, (3. Furthermore, we have the following. 
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Lemma 3.8 For any a G Z>q and positive integer N , the function lN,a is a polynomial in 
N of degree [||q;||/2] with coefficients in C{P) fl C°°{P°). More precisely, lN,a is of the form 



[||"ll/2] 
1=0 



(41) 



where pa,i G C{P) nC°°(P°) and Pa,i does not depend on N. Furthermore, D^pa^i E C{P) fl 
C°°{P°) for every (3. In particular, we have lN,a{^) ~ 0(A^["""/^1) uniformly on P. 

Lemma [3l8] can easily be shown by induction on k = \\a\\ with Lemma [321 Note that the 
functions Pa,i can be computed inductively. Indeed, by substituting (HTi) to (HUj) . we have 



ll"ll+l)/2] 



(=0 



DjP^^l{x) + ^aiIi^e,+e^{x)Pa-ea-l{^) 



i=l 



(42) 



3.2 Proof of Theorem [2791 



By using the Taylor expansion of / around x E P, the function B^{f) is represented as (13T]) . 
(!32l) . The error term S2n,N in (1311) is estimated as (l33l) . and hence, by Lemma [378| we have 
II 5'2„,Ar II c(P) = 0(||/||c2n(p)A^""). Thus we have 



a <2n-l 



r/7V,a(x)+0(||/||c-(p)iV-"). 



Then, by substituting (HTl) for the above, we have 

1^ 



|ja||<2n-l 



n-1 / 2v 



a! 

2n-l / 2n-l 



N~ 



v=n \ 1=1/ \\a\\=l 



n-1 / 2v 

E EE 

v=Q \ l=u \\a\\=l 



/W(x) 



Po.,i-u{x) iV-'^ + 0(||/||c--i(P)iV- 



Therefore, if we set 



2u ^ 



(43) 



l=u ||a||=; 



we obtain 



n-1 



BAf){x) = J2N-''Ux,d)f + 0{\\f\\c^.^P^N~^), 



i/=0 
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which shows the asymptotic expansion (pTI) . The formulas for the differential operators 
Lo(x, 9) and Li{x,d) are obtained from Lemma [3171 and (143|1 . 

Next, we show that the expansion (pT!) can be differentiated any number of times. First of 
all we note that, if djB^{f) has an asymptotic expansion of the form ^^>q N~^gjy{x) locally 
uniformly in x £ P° for each j = 1, . . . , m, then it is easy to show that gj^^ = djLi,{x, d)f. 
Instead of using the partial differential operators dj, we use the operators Dj defined in ( 15^ . 
From Lemma 13.21 we have 

D^B^if) = N{Br,{X.jf) - X.B^if)), that is, D.B^ = N[Bn,X^]. (44) 

Therefore, since Lo(x,9) = 1, we have DjB^^f) ~ Yliuxd^^^i^u+ii^id), Xj\f . This holds 
uniformly on P. From this we also have DjL^{x,d) = [L^^i{x,d),Xj]. Next, assume that 
we have 

D'^BM) ~ E N-'D-L^ix, d)f (45) 

uniformly on P for any / G C°°{P) and a G Z>q with ||a;|| < k where is a fixed positive 
integer. Then, by (Hlj), we have 

D.D^BnU) = D'^D.BnU) = ND^{BN{X,f) - X^B^U)). 

and hence 

D^D^BnU) ~ 5^iVi-'^(/^"L,(x,a)X, - D''X,L,{x,d))f 

= J2N-'D'^[Lu+iix,d),X,]f = J2N-''D'^D,L,ix,d)f. 

Thus, by induction, the asymptotic expansion (H5|) holds for any a uniformly on P. Since 
A[x) is invertible on P°, we conclude the assertion. □ 



Remark 3.9 One may think that the theorem can be proved by using the method of sta- 
tionary phase. Indeed, we have the formula 

BN{f){x) = / mm^iO =(^Y [ e-^*(^'^'«)/(?/) (46) 

where / is the Fourier transform of an extension of / G C'^{P) to the whole space M"* as a 
compactly supported smooth function which is also denoted by /, and the function ipN,x is 
the characteristic function, 

of the probability measure dB^ . Since dB^ is defined by the iV-th convolution power of the 
Bernstein measure i3, we have 

VnAO = ^{^,aNf, ^{x,0 ■■= J^e^'^dB^iz). 
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Then, the function y,^) in the right hand side of (146|) is given by 

^{x, = ^( O - log vi^^ 0- 

Thus, one would find asymptotic expansion for B]^{f) by using the method of stationary 
phase (for example, using Theorem 7.7.5 in |Holj ) with the phase function $(a;,?/,.^). Note 
here that, in this computation using stationary phase method, one might not need to assume 
the condition (2) in Definition 12.11 for the section B : P —>■ M.{P). Theorem 7.7.5 in [Hoi] 
gives us an effective formula for each term of the expansion. However, it is not still quite 
easy to compute each term of the expansion. Instead of using the method of stationary 
phase, we used in the above proof a beautiful idea given by Hormander ( [Ho2j ) which gives 
a computable representation fl^T]) . fH21) . fH5|) of each differential operator L^{x,d) when B is 
a Bernstein measure. 



3.3 Proof of Theorem [230] 

We keep the notation described before the statement of Theorem 12.101 in Section [2l First, 
we note that the functions ms,c,a, a & S defined in (j26l) are continuous up to the boundary 
dP of P. Before giving a proof of this fact, we describe the values of ms,c,a at a; G dP. Let 
K he a (relatively open) face of P. Let 

Xk = spanjgja — P ; a, P & S H K}, 

and let C M"^ be the annihilator of Xk- Then, define the map hk : X —>■ K hj 



(47) 



If t,t' E X satisfy t — t' E X^, then /ix('r) = /i/^(r'), and hence the above defines a map 
Hk '■ Xk = W^/X^ — * K, and it is a diffeomorphism from Xk onto K ([F]). For any 
u e M"*, we define \{u) = minj^gp( y, u ). We fix u G X^ satisfying 

K = {yeP; {y,u) = X{u)}. 

Then, for any x = /i/<(T) G K with r G R™, we have lim^^+oo f^s,c{T — tu) = x and 

ms,cAx) = 1™ "^5,c,a(Ai5,c(r - tw)) 

t — ^+oo 



a^SnK, (4^ 

a e SnK. 



Lemma 3.10 The functions Tns^c,a{x) (a G S) are continuous on P. 

Proof. Let K C dP be a relatively open face of P. Let Xn G P° and x G K 
satisfy — a; as n — > oo. We show that ms^c,a{xn) — > rns^c,a{x). (For the case where 
Xn is contained in a face L ^ P°, one can discuss as in the following with replacing P° 
by L.) Let := rs,c(a;„) G M™. Take u e X^ such that F = {y G P ; ( y, ?i ) = A(m)} 
with A(m) = minygp( ?/, u ). For any A C 5*, we define Xa(t) = J2a(^A'^('^)^^°'''^^ ■ 
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Sk := S n K and ck = niinQ,g5\5^{( a, u ) — \{u)}. We note that {a,u) — \{u) = if and 
only if a G Sk, and hence ck > 0. Since < ( x„, u ) — X{u) — > as n — >• cx), we have 

< CK /ff ^"""^ , , < ( fisArn), u ) - X{u) ^ 0, 



and which shows 



0, — — — r > \) [n ^ OO). (49) 



XSk i^n) + XS\Sk (^n) XSk i^n] 

From this, we have hm„,_,oo fns,c,a{xn) = = ms,c,a{x) when a G S\Sk- Decompose r„ G M"^ 
according to the decomposition = Xk © -^^^ as = + ^n, G X/^, Un G X^. Then, 
I9l) also shows 

Ec(a)e^°''^"^ 
^ — —a^x (n oo). 

Since /ik '■ Xk ^ K is a diffeomorphism, we take ^ G Xk such that fiKiO — ^- Then, 
the above means lim^^oo /^_ft:(^n) = /^^^(O) hence we have lim„^oo^n = ^- From this we 
have, for a E Sk, 

. . c(a)e<"'«"> ^ c(a)e<"'«> 
as 77, ^ cxD, which shows the assertion. □ 



Remark 3.11 For the original Bernstein polynomials B]\f{f) defined in ([T]), one has B]\f{f){0) = 
/(O), i?Ar(/)(l) = /(I). Our Bernstein approximations B^if) defined by a finitely supported 
Bernstein measure Bs,c also have similar property. In fact, by fj48|) . if x is a vertex of P, 
one has ms,c,x{x) = 1, ms,c,a{x) = when a ^ x. This combined with f|T2l) shows that 
BN{f){x) = f{x) when x is a vertex of P. 

Now, we define the function 6s,c G C°°{P°) by 

55,c(x) =log (5^c(«)e<°'^«-(^)M - (x,rs,e(x)), (50) 

\aeS / 

where, as in Section [2], ts,c '■ P° —>■ is the inverse map of the diffeomorphism fis,c '■ 
P°. Then, we have 

ms,cA^) = c(a)e-^s,c{x)+(a-x,rs,.{x)>_ ^5^^, 

Completion of proof of Theorem I2.1UI As in |TZj . we have V5s,c{x) = -ts^x) 
and V^Ss,c{x) = —As^x)^^ for x G P°, where As^c is defined by 

AsA^) = '^'m'S,cA^)i'^ - ^) ® (oi-x), (52) 
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which is non-degenerate on P° (see |TZ] ). Then, it is not hard to show directly that the 
measure (l28l) defines a Bernstein measure on P with the defining matrix As,c{x)~^- 

Next, let us prove the converse. Let B : P ^ A4{P) he a finitely supported Bernstein 
measure, and let S = supp(i3(x)) for some, and any x G P°. Let K : P° ^ Sym(m, R) be 
the defining matrix of B. We write B{x) = X^^^g ?tIq,(x)5q,, x G P. Then, by Lemma [XT| 
the matrix K{x) is non-degenerate. Now, fix an arbitrary b* G P° and > 0. Define the 
function 6b* G C°°(P°) and the map n* : P° ^ by 

n*{x)= f K{b* + six-b*))ix-b*)ds, 
Jo 

5b*{x) = \ogV- [ {n*{b* + s{x-b*)),x-b*)ds. 
Jo 

Lemma [3751 shows that {dTb*)x = K{x), V5b*{x) = —Tb*{x), x G P". Since ma{x) > on P°, 
we define ga{x) = log ma (x) — {a — x, Tb*(x) )• Then, since Vga = U*, there exists a constant 
c{a) > such that ga = —5b* + logc(a) for each a & S and hence 

m«(x) = c(«)e-^''*(^)+<"~^'"''*(^)\ (53) 

Taking the sum of ( l53i) over a G S", we have 

5b* {x) = log ( 5^ c(«)e< M - ( x, n* (x) ) . 

Since Tb*{b*) = 0, we have rriaib*) = c{a)/V. It follows from this and the identity Xlaes' — 
b* that ^ 

^ = ^c(a), b* = -^c{a)a. 

From these formulas, we obtain 6b* (b*) = logV. It is not hard to show that the map Tb* is 
injective and hence is a diffeomorphism from P° to the image of Tb* . (Note that the differential 
of Tb* is K, which is positive definite on P°.) Denote its inverse by /x;,* : Im(rf,.) — >■ P°. 
For any x G P° and r G Im(rfe*), we define /r(x) = (x, r) + Sb*{x). Then, we have 
V/r(x) = T — Tb*{x) and V^/t(x) = —K{x). Thus, the point x = iib*{,T) is a unique critical 
point of /r, and /r attains its maximum there. Then, we obtain 

frix) < f^{fIb*{T)) = logX5,c(r), 

where the function xs,c, defined by 

X5,c(r) = 5^c(a)e<---\ (54) 

aeS 

depends only on the constants c(a) > and the set S. In the above inequality, the equality 
holds if and only if x = fib*{T). Therefore, we obtain 

6b* (x) < log xsAt) - 
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for every r G Im(rf,.) and equality holds if and only if r = Tb*{x). Then, as in |TZ] . we 
obtain 

5sA^) = inin (logX5,c(i") - {x,r)) < 5b*ix). 

Since Im(r5.) C is open, the point r = Tf,*{x) is a local minimum of the function 
T log Xs',c(t) — ( X, r ) on R™. Since, the local minimum of this function is unique, and it 
is given by r = ts^x), we conclude r^. (x) = ts^x) for every x G P°. This completes the 
proof. □ 

Remark 3.12 In the proof of the converse direction of Theorem 12.101 above, there are no 
restriction for the weight function c : S* — >■ M>o, which is defined as an integral constant. 
This is because any choice of the weight function c : S ^ M>o can define a Bernstein 
measure. However, two different weight functions might give the same Bernstein measure. 
In fact, when P is a standard simplex and S = P (1 Z™, the set of vertices of the simplex 
P, any choice of the weight function produce the Bernstein measure given in Example 12.61 
In general, let A C P be a finite set whose convex hull is P. Let C Ai{P) denote the 
convex hull of the Dirac measures 6a with a & A. Then, the image of the finitely supported 
Bernstein measures B : P ^ Ai{P) with the support S are contained in A^, and the image 
of the dilated convolution powers dB!^ are contained in A±g . When P is a standard m- 



dimensional simplex and S = P H Z™, As is also a simplex with the same dimension, and 
the restriction of the barycenter map to A^ is a diffeomorphism between As and P. Hence, 
in this case, there are only one section of the barycenter map whose support (in the sense 
of Definition 12.31) is S. 



3.4 Large deviations principle 

In this subsection, we discuss the large deviations principle for finitely supported Bernstein 
measures. Throughout this subsection, let B = Bs^c '■ P —* A^(P) denote the finitely 
supported Bernstein measure with the support S G P and the weight c : S ^ R>o- For each 
a E S, the coefficient of B is denoted by ms^c,a{x) defined in 
For any a; G P, we define a function Xx = Xs,c,x by 



T G W^. (55) 



By using the function xs,c defined in ( 154|) . we find that 



= ^^^'-'-f (56) 

Lemma 3.13 We fix x E P . Then the dilated convolution powers dB^ defined in ([8]) from 
the finitely supported Bernstein measure B = Bs,c satisfy the large deviations principle with 
the speed N and the good rate function given by 

P{y)= snp{{y,T)-A^{T)}, A-(r) = logx.(r). (57) 
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Proof. Consider the infinite product of dB^ on tlie infinite product space f2 = P x 
• • • X P X ■ ■ ■ . Let Xj^x : f2 — >■ P denote the projection in the j-th coordinate. Then Xj^^s 
form a sequence of independent identically distributed random vectors with the distribution 



dBx- Then, the distribution of the empirical means S^^^ '■= ^J2f=i-^j,x is given by the 
dilated convolution powers dB^ . Thus the assertion follows from Cramer's theorem ( |DZj ). 
□ 

Next, we describe the rate function (|57j) by using (H71) . ( HHl) . To describe the rate function 
(l57j) more concretely, we need the fact that the function 6s,c in ( l50l) is continuous on P. Let 
K he a (relatively open) face of the polytope P. Let tk '■ K ^ Xk be the inverse of the map 
liK '■ Xk — > K defined in f H7|) . When K = P°, we have fipo = fis,c and Tpo = ts,c- As in the 
proof of Lemma l3.10\ for each A C S, we set Xa(t) = J2aeA'^('^)^^°''^^ = S H K. 

Define a function 6k on K by 

6Kix) = \ogxsKiTKix)) - {x,TKix)), X E K. (58) 

Note that, we have 6po = 6s,c- 

Lemma 3.14 The function 6s,c is continuous on P, and its restriction to each face K is 
given by 6k- 

Proof. As in the proof of Lemma [3.101 take u e X^ such that K = {y E P ] {y,u) = 
X{u)} with X{u) = mmz<zp{z,u). Then, first of all, we claim that, for x E K, the following 
holds: 

lim 6s,cilJ'sA^Kix) - tu)) = 6k{x). (59) 
Let us prove ( l59|) . It is easy to show that, for any x E K, 

\{fisArK{x)-tu)-x),u\ < Ce-'^K ^s\sArK{x)) ^^^^ 

XSk{^k{x)) 

where we set Ck '■= nain^esys^^, ( a, u ) — A(m). For any ^ G Xk, we set 

< P^(t,0 := Yl c(a)e<"'^>-*«"'">-^(")) < e-'^-xs\sAO- (61) 

Since x E K, we have {x,u) = \{u), and hence 

SsAf^sA'^Kix) - tu)) 
= \ogxsKi'TKix)) - {hsA^k{x) -tu),TKix))+ (62) 

+ t{lJ,sA'TKix) - tu) -X,U) + \og{l + RK{t,TKix))/xSKi^Kix))). 

Now, by (l60l) . (16T|) and the fact that UsA'^k^x) — tu) tends to x G -ft' as t ^ +00, the right 
hand side of fl62l) converges to 6k{x), which shows ( 1591) . 

Next, we take x E K, and {x„} C P° such that Xn ^ x as n ^ 00. Let ^ = rx(x) 
and r„ = ^„ + M„ = TsAxn) with ^„ G X^, £ In the proof of Lemma [3.101 we 

have proved that ^„ — » ^. Let PK^n) = {z,Un) with z E K, which does not depend on the 
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choice oi z G K. We note that X5k(^" + — ^^^^^^XsKi^n)- From this and (H9!) . for each 
a & S \ Sk, we have 

c(a)e<"'5n>+(a,«0-PK(n) < e-^-(")xs\SK(a + Wn) 0, 

and hence {a,Un) — Pxin) — oo as ^ oo for each a G 5 \ 5*/^. Now, by using this fact, 
we have 



logX5K(^n)- 2^ -—{{a,Un)-pKin)) + o{l)-^\ogxsKiO- 



a£S\S_ 



In the above, one can compute the term {xn,Un) by using the relation x„ = fJ^sA^n + 
and the definition (125|) of the map /zg^c- Since 5s,c(^n) = —{xn,^n + Un) + log x^.tl^n), we 
conclude that Ss,c{xn) — > In case where is contained in a face L 7^ P°, one can 

discuss in the same way as above with replacing P° by L to conclude the assertion. □ 



Remark 3.15 The function 6k on the face K defined in f lSS]) is continuous on K. This can 
be shown in the same way as in the above proof. The restriction of Sk on a (relatively open) 
face L of ii" is given by 5^. 

Proposition 3.16 Let K be a relatively open face of P. Letx G K. Then, the rate function 
P{y) is given by the following: 

I SKix) - 5K{y) + ( X - tk{x) ) yeK. ^^""^ 

Proof. Let /sT be a face, and we fix a; G /sT. By (iH]), (ESD, we have Xxir) = ^''^^T'^tr^ - 
Thus, we have 

I^'iy) = JK{y) + 5K{x) + {x-y,TK{x)), 

My)= snp{{y,T)-logxsAr)}- ^^^^ 

This holds for any y G M™. For any y G M™, we set 

Ixiy) = snp {{y,^) - hgxsKiO}- 

In (16^ . we decompose t = ^ + u with ^ G X^, u G X^. Then, for y E K, we have 
+ -logXx(^ + M) = {y,^) -\ogXxiO and hence = JKiy) for each ?/ G iT. 

Now, first of all, we assume y E K. Then, for ^ G Xk, we have logX5K(0 = f^K^O^ which 
shows that ^ = TK^y) is a unique critical point of the function ^ i— >• (y,^) — logx^^lO- The 
Hessian of the function log xskIO is non-negative, and hence we have 

Ixiy) = {y,TK{y)) -hgxsKi'^Kiy)) = (65) 

which shows ( l63l) for y E K. We set F^{y) = 5^(3;) — 5x(z/) + {x — y,TK{x)) for y E K, which 
is continuous in ?/ G -ft'. Since /^(y) is lower semi-continuous in y, we have F^{y) > P{y) for 
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any y E K. Now, assume that y is contained in a face L of K, and we show the inequahty 
Ffciy) < I^iy) foi' such a point y, which imphes (!63|) . Since /^d/) = Jxiy) for any ?/ G it 
suffices to show that Ixiy) > —h<{y)- Take m G X^;^ such that L = {z E P ] {z,u) = \{u)} 
with A('u) = min2gp( z,u). Then, for any C, G and t > 0, we have 

»&Sk\Sl 

Note that Rxit, ^, m) — > as t — ^ +cxd. Since y E L, we have A(u) = {y,u), and hence 

where 0^(1) denotes a term tending to zero as t — +oo. Since Ixiy) = Jxiy) > {y,^ ~ 
tu ) — logX5x(^ ~ ^'") ^ ^ -^L and t > 0, we obtain 

/i.(i/)>(i/,0-iogx5,(0, eeXz.. 

Since the supremum over ^ G Xj;^ of the right hand side above is = —dxiy), we 

conclude Ixiy) ^ ~^K{y) for y E L. Finally, we show that I^{y) = +oo when y ^ K. We 
have 

(y,^) -logXSxl^) = -logX5K(0) + / {y - l^K{tr),T) dt. 

Jo 

Note, in the above, that /i/^(tr) G -ft' for any t E [0, 1] and r G M™. Since y ^ K, we can 
take M G such that {y — z,u) > for any z E K. Then, by setting t = ru with r > in 
the above, we conclude that Jxiy) = +oo and hence I^{y) = +oo. □ 

Remark 3.17 When the polytope P is a lattice polytope satisfying Delzant condition, the 
formula for the rate function fl63p coincides with the formula given in |SoZ2j . Proposition 
5.2 for the rate function in a large deviations principle of the Bergman-Bernstein measure 
(see Section Hj) defined by the Fubini-Study metric on the toric manifold obtained through 
a monomial embedding. Indeed, for example, when x E P°, the rate function P{y) is given 
by 

i^'iy) = SsA^) - ^sAv) + i^-y, ^sA^) ) 

= -SsAv) + ^ogxsA^sA^)) - iy^'^sA^) )• 

Then, the function log X5_c is a Kahler potential for the Fubini-Study metric and —Ss,c is its 
Legendre dual. See [Z], |SoZ2j and the following section for details. 

4 Bergman-Bernstein approximations 

In [Z], Zelditch introduced the notion of the Bergman-Bernstein approximations for functions 
on Delzant polytopes P (a lattice polytope with the property that each vertex of P has 
exactly m = dim P edges and m lattice vectors incident from the vertex along the edges form 
a Z-basis of the lattice). We explain here this notion for projective toric manifolds obtained 
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by monomial embeddings and the difference between Bergman-Bernstein approximations 
and Bernstein measures defined in this paper. 

Throughout this section, we assume that the polytope P is Delzant with vertices in 
Z"". We set S = PC] IT^ and = NP ^ Z™. We fix a function c : S ^ R>o. Let 
^s,c '■ (C*)™" — > CP'"^'^^ be the monomial embedding defined by 

Then, the toric variety Mg c is defined by the Zariski closure of the image of the monomial 
embedding $5,^ : (C*)*" ^ CPl^l^^ f [GKZ] l Note that, in general, Ms,c may have singular- 
ities. However, it is well-known ( [GKZ] ) that, if P satisfies the Delzant condition, then M^c 
is a non-singular compact Kahler manifold. Consider the (restriction of the) Fubini-Study 
Kahler form ups on Ms^c- Then, the action of the real torus T"^ C (C*)"^ is Hamiltonian with 
respect to the symplectic form cups- On the open orbit (C*)"^ = $5^c((C*)™'), one can take a 
T™- invariant Kahler potential (p. Since is a function on (C*)™ invariant under T'^-action, 
it defines a function on M™, which we denote by 'ps,c- Let /i : Ms,c — > P be the moment map 
of the Hamiltonian action of on {Ms^ci'^fs)- Since /i is also T'"-invariant, it defines a 
map ixs,c '■ — *■ P, which is a diffeomorphism onto the interior, P°, of P. 

The symplectic potential associated to ip is the Legendre dual of the function ips,c 
associated to the Kahler potential which is defined by 

U^X) = ( X, Ts,c{x) ) - VsAtsAx)), X e P°, 

where ts^c '■ P° is the inverse of ns,c- The Bergman-Bernstein approximation z/7v(/) 

(the notation BhN{f) is used in [Z], where h denotes a Hermitian metric on the hypersection 
bundle Lsc '■= 0{1)\ms over Msc whose curvature is o^ps) of a function / on P is defined 
by ' 

1 ^ N{u^{x)+{'y/N-x,Ts,c{^))) 

Mf){x) = - fil/N) -—^ , (66) 

where Qh'^il) is the squared L^-norm of the monomial with weight 7, which is regarded 
as an element of H%Ms,c,L^^), the function H^ {z, z) is the Bergman-Szego kernel for 
H"{Ms,c,Lll^), and z e (C*)™ satisfies /i(z) = x. 

To compare the Bergman-Bernstein approximation fl66|l with our Bernstein measures, we 
take the Kahler potential 

cp{z)=\ogY,cia)\z-\\ ze{C*r. 

of ujps on (C*)™". Then, the corresponding function (ps^c coincides with the function xs,c on 
M"* defined in (IMI) . Therefore, the symplectic potential coincides with the function —Ss,c 
defined in (ISUl) . The quantity Qh^ipf) is, as in [Z], given by 

Note that the restriction of the moment map /i to the open orbit is given by 
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Thus, the map j^s^c '■ 1^™ P° induced by the moment map /i is nothing but the map 
defined in fl25l) . 

Lemma 4.1 Let m^^(x) he the function on P defined by (fT2l) {with rUa replaced by ms^c,a 
defined in (!26l) ). Then, the Bergman- Bernstein approximation (!66l) is written as 

where the quantity -Rat (7) anc? t/ie function Ilfyf{x) are given by 

Rwil) = J ml Nix) dx, 

Proof. First, we note that the function ms^c,a{x) is written as flHTl) . Thus, the functions 
ml j^{x) defined in (fT2|) is written as 

where "Pat (7) = Vs^cnIi) is the weighted number of lattice path, 

^7v(7):= 5^ c(/3i)---c(/5;v). 

/3i,...,/3jvg5;/3i+-+/3jv=7 



From this, we have 
which shows that 



RN{l)=QhHl)'PN{l), (69) 



Qh^ii) Rn{i)' 

Note that the Bergman-Szego kernel IIn{z,z) with z G (C*)™ is written as 

n^(^, ^) = E oM^, = n^(^)' (70) 

where x = fi{z), and hence we have the assertion. □ 

Therefore, it would be natural to call the probability measure 

■= E ^^irr^w^ (71) 

the Bergman-Bernstein measure. (The measure duf^ defined above equals the measure 
with ^{z) = X in [Z].) Note that the Ber gman-Bernstein measure duff is not a section of the 
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barycenter map b : Ai{P) P. In fact, by Lemma [3^ and Theorem 12. 10[ it is easy to show 
the following formula: 

j^DlogUN{x) = b{du-^)-x, (72) 

where, for / G C°°{P) and a; G P, we set Df{x) := As,c{x)V f (x) with the matrix As,c{x) 
defined in ( l52l) . and b{di'%) denotes the barycenter of the Bergman-Bernstein measure dufj. 
Then, it is natural to ask when the Bergman-Bernstein measure duf^ coincides with the 
dilated convolution power dB^ of the Bernstein measure Es,c{,x) defined by (1281) . For this 
question, we have the following proposition. 

Proposition 4.2 Let dv^ he the Bergman-Bernstein measure defined by (17T|) . and let dB^ 
denote the dilated convolution power ([8]) of the finitely supported Bernstein measure Bs,c{x) 
defined in (!28l) . Then, the following four conditions are equivalent. 

(1) Rn{i) is constant as a function on the finite set S^. 

(2) The function IIn{x) defined in fl68|) is constant. 

(3) The barycenter of duff is x for each x E P. 

(4) duff = dB^ for each x E P. 

Remark 4.3 By fl70|l . we know that the function IIn{x) on P is the function induced by 
the restriction of the Bergman-Szego kernel to the diagonal. Hence, the condition (2) in 
Proposition 14.21 is equivalent to that the Fubini-Study Hermitian metric on i^f ^ is a 
balanced metric ([Dj). 

Proof. According to the formula ( 1721) . it is obvious that the conditions (2) and (3) 
are equivalent. Assume that the condition (1) holds. Since dB!^ = ^^^Sn ™'s,Ni^)^'y/N' is a 
probability measure. Hat is constant because of its definition (!68|) . which shows (2). Next, 
assume that the condition (4) holds. Then, we have Hat (x) -Rat (7) = 1 for each 7 G Sn and 
X E P. Since Un does not depend on 7, the condition (1) holds. Finally, assume that the 

condition (2) holds. We set mj^^(x) = njv^fliv(7) ^^^^ = Yl,'yi^SN'^'s,N{^)^i/N ■ Then, 
by Lemma [3. 2 [ we have 

Vml^{x) = ml^{x)K{x){^ - iVx), (73) 

where K{x) = As^c{x)~^ is the defining matrix of the Bernstein measure Bs^c '■ P M.{P). 
Let T/v : C{P) C{P) be the linear map defined by duf^, that is, 

T;v(/)(x) := rhl^ix)fil/N), f E C{P). 

Then, by (173|) . we have 

D.T^if) = NT^iX.f) - NX.T^if), f E C{P), 

where Xj is the multiphcation operator, {Xjf){x) = Xjf{x), and Djf{x) is the j-th compo- 
nent of As^c{x)V f{x). As in the proof of Proposition 13. 6^ define the first-order differential 
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operators Ltvj (j = 1, . . . , m) by L^^j = Dj + NXj. Since [Lnj, L^^k] = for each j, k, we 
write = ■ ■ ■ i^^'"^ for a = (ai, . . . , am) G Z>q. Then, as in the proof of Proposition 
13. 6[ we have 

L%Tn{xp) = iVll"llTjv(Xa+/3), and hence T^(xa) = " 1, 

where Xa denotes the monomial with weight a. Note that the last expression in the above 
depends only on the matrix As^d^), and the same formula holds for the Bernstein approx- 
imation B]\r{xa) instead of T/v(Xo)- Therefore, we have T^{f) = B^i^f) for any / G C{P), 
which shows the condition (4). □ 

One of advantages of using the Bergman- Bernstein approximation i^7v(/) is that it satisfies 
the identity 

f nj,{x)u^{f){x)dx= Yl /(^/^)- (74) 

Then, once one find an asymptotic expansion of z^Ar(/) as ^ oo, one would have an 
asymptotic expansion of the Riemann sum by using a well-known asymptotic behavior 
of the Bergman-Szego kernel function Il]sr{x). This is the idea used in [Z]. (More precisely, 
Zelditch obtains an asymptotic expansion of the numerator in (l67|l .) For the Bernstein 
approximation B^lf) introduced in this paper, we can not, in general, expect that the 
identity like (17i|) holds for B^if)- Instead, our Bernstein measures can be used for general 
polytopes. 

It might be possible to find asymptotic behavior of the Bergman-Bernstein measure 
defined by fl67|) even for general polytope P. However, for this, one might need to analyze 
in detail the behavior of the functions rug ^{x) when and x are close to the boundary 
of the polytope. 
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